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THE HANNAY-BERRY PHASE OF THE VIBRATING RING 
GYROSCOPE * 


SEAN B. ANDERSSON t AND P.S. KRISHNAPRASAD * 


Abstract. 

In an analysis published in 1890 G.H. Bryan investigated the retrograde precession of the nodal 
points in a vibrating, rotating shell and wrote down a formula relating the rate of precession to the 
rate of rotation. This effect has been utilized in the design of various vibratory gyroscopes includ¬ 
ing modern MEMS-based devices. Existing analyses model these systems with a pair of harmonic 
oscillators coupled through the Coriolis force (the normal mode method). In this work we utilize 
the theory of moving systems developed by Marsden, Montgomery, and Ratiu to show that the 
nodal precession can be understood as a geometric phase with respect to the Cartan-Hannay-Berry 
connection. This approach allows us to explicitly characterize the simplifications of the linearizing 
assumptions common to previous analyses. Our results match those of Bryan for small amplitude 
vibrations of the ring. We use the inherently nonlinear nature of the moving systems approach to 
calculate a (small) correction to the rate of precession of the nodes. 

Key words, geometric phase, holonomy, fiber bundle, vibratory gyroscope, MEMS, rotation 
sensing 

AMS subject classifications. 53C29, 70G45, 70H05, 70K65,74F99, 74K10, 81Q70, 93C95 

1. Introduction . In 1890 G.H. Bryan published a paper on the nature of the 
beats generated when a vibrating shell is rotated about its central axis [8]. The 
phenomenon he describes is quite easy to observe; simply take a wine glass, strike it 
to produce a clear tone, and then rotate it about its stem to produce audible beats. 
Bryan noticed that these beats are the result of a precession of the nodal points with 
respect to the shell itself and provided the following reasoning. Consider a ring or 
cylinder rotating about its central axis and vibrating with nodes at B,D,F, and H as 
indicated in the left-side image of Figure 1.1. The material points at A and E are 
moving towards the center O. This increases their actual angular velocity above that 
of the imposed rotation and gives them a relative angular acceleration in the direction 
of rotation as represented by the arrows at A and E in the right side image of Figure 
1.1. Similarly, the material points at C and G are moving outwards and thus their 
angular velocity is reduced. Those at B and F are moving with greater total angular 
velocity than the rest and thus experience a relative outwards acceleration due to a 
greater centrifugal force. Finally, the material points at D and H are moving with the 
least angular velocity and thus experience a relative acceleration inwards. Comparing 
the arrows in the two images of Figure 1.1 reveals that the effect of these relative 
accelerations is to cause retrogade motion of the nodes relative to the ring. Using 
classical variational techniques Bryan derived a linearized partial differential equation 
describing the behavior of the system, found a formula for the rate of precession and 
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Fig. 1.1. Nodal precession in cylinder or ring (Figure from [8]) 


discovered that this rate is proportional to the rate of platform rotation. 

Due to the immense number of potential applications, research on gyroscopes has 
been active for many years. Devices have been proposed, analyzed, and produced 
using a variety of materials and techniques. Because they lack rapidly spinning parts, 
have low power requirements, and are inherently scalable, vibratory gyroscopes have 
become particularly popular [21]. One of the most successful initial designs was 
Delco’s Hemispherical Resonator Gyroscope (HRG) [15] due to Loper and Lynch, 
which was able to achieve performance levels equal to the best ring laser gyroscopes. 
This design, shown diagrammatically in Figure 1.2, consists of a quartz hemispherical 
resonator supported on a central stem and contained inside an evacuated housing. 
As predicted by Bryan, the nodal points of the vibration in the hemisphere precess 



Fig. 1.2. Delco Hemispherical Resonator Gyroscope (Image from [14]) 


with respect to the shell as the device is rotated. The HRG is driven into elliptical 
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vibration as shown in Figure 1.3 and the resulting precession rate is about 0.3 of 
that of the imposed rotation. The device operates over a wide temperature range, has 
high operating acceleration ranges, low accleration sensitivies, and negligible magnetic 
sensitivity [16]. Similar ideas were used in the design of a vibrating disc gyroscope by 



Fig. 1.3. Precession of Delco HRG (Image from [14] following [15]) 

Burdess and Wren [10]. 

With the explosion of MEMS technology constant innovations are resulting in 
smaller, cheaper, and more accurate devices. Existing MEMS-based devices include 
tuning-fork [5] and vibrating-ring designs [4, 22] such as the one shown in Figure 
1.4 (provided by Douglas Sparks of Delco Automotive Systems). Additional designs 
proposed include a vibrating cylinder [9, 28] and a surface acoustic wave generator 
[27]. 



Fig. 1.4. MEMS Gyroscope (Image courtesy of Douglas Sparks of Delco Automotive Systems) 
These gyroscopes all take advantage of the same physical effect, the Coriolis force 
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arising from the non-inertial character of the rotating frame of the system [14]. Mod¬ 
ern analyses are linear in nature and view the Coriolis force as providing a coupling 
between two vibratory modes of the structure [25]. It is desirable, however, to have a 
method which, at least in principle, can be extended to a nonlinear theory and which 
provides a unified setting for understanding a variety of systems in which the Cori¬ 
olis force plays a role. Motivated by these considerations we are led to an approach 
developed by Marsden, Montgomery, and Ratiu based on modern developments in ge¬ 
ometric mechanics. This method is known as the moving systems approach [18]. The 
technique, descending from the classical work of E. Cartan [11], describes the effect of 
imposed motion on a system as a geometric phase with respect to a particular Ehres- 
mann connection called the Cartan-Hannay-Berry connection. This geometric phase 
is termed the Hannay-Berry phase. In previous work we have used this approach 
to determine that the Hannay-Berry phase in a rotating, equal-sided, spring-jointed, 
four-bar mechanism is zero [3]. (See also [2].) 

Interest in the effect of geometric phases in physical systems was spurred by work 
of Michael Berry in the early 1980’s [6, 7]. Berry’s initial result, surprising in its 
simplicity, has sparked a great deal of research that continues to this day (see, e.g. 
the compendium volume of Wilczek and Shapere [24]). 

While existing techniques have proved to be effective, as evidenced by the various 
devices constructed from those principles, it is to be expected that a deeper under¬ 
standing will emerge by appealing to a nonlinear, geometric approach directed at more 
accurate constitutive models. The fact that Bryan’s result is geometrical in nature 
is intriguing and suggests that these modern tools may prove useful. In this paper 
we begin this study by applying the moving systems approach to the vibrating ring 
and show that Bryan’s result can be interpreted as the geometric phase. The main 
contribution of this work is to provide, through the use of geometric techniques, a 
clear understanding of the role of the linearization assumptions common to previous 
analyses and to make explicit the simplifications that result. Furthermore we calcu¬ 
late the effect of the neglected terms by deriving a formula for a correction term to 
the rate of precession of the nodes of vibration. Finally, this work illustrates, through 
the reasonably complex example of the ring gyroscope, how the methods of modern 
geometric mechanics are useful in understanding physical effects in systems which are 
useful for sensing and as a result hopefully motivates new applications of the ideas. 

The remainder of the paper is organized as follows. In the next section we provide 
a brief background the proper mathematical setting for geometric phases, namely 
that of connections on fiber bundles. In Section 3 we summarize the moving systems 
approach of Marsden, Montgomery, and Ratiu. In Section 4 we apply the technique 
to the vibrating ring and derive the Hannay-Berry phase which is then compared in 
Section 5 to the results of Bryan. The imposed rotation on the ring gives rise to 
nonlinear terms in the equations of motion even when a linear approximation of the 
nominal (non-moving) ring is used and in Section 6 we calculate a correction term 
incorporating these effects. We conclude in Section 7 with a few comments on the 
results and on future work. 

2. Background . The natural mathematical framework from which to approach 
geometric phases is that of connections on fiber bundles. In this section we present 
a review of the necessary background on fiber bundles, connections, and holonomy. 
Additional references for this material include [12, 19] and references therein. 

2.1. Fiber bundles and Ehresmann connections. A fiber bundle is de¬ 
fined as follows. Let P, F, and B be manifolds referred to as the total space (or 
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bundle space), the fiber, and the base space respectively. Let 7r : P —> B be a 
surjective submersion. We require that P be locally a product space, that is, for every 
b £ B there is a neighborhood U of b such that 7r _1 (f7) is diffeomorphic to F x U. 
The fiber over b, 7r -1 (i>), is a diffeomorphic copy of the fiber F for every b £ B. See 
Figure 2.1. The bundle is denoted by the triple ( P,F,B ) or by the projection map 



a 


Fig. 2.1. Fiber bundle, fiber over b, and splitting of the tangent space 

7 t : P ^ B. If the bundle is globally a product bundle, P = F x B, then it is called a 
trivial fiber bundle. 

Given p £ P, there is a natural subspace of T p P (the tangent space to P at p) 
called the vertical space at p, denoted by V p and defined by V p = kerXpTr. Here T p n is 
the linearization of the projection map evaluated at p. The union of these subspaces 
over all p is called the vertical subbundle V, i.e. V = U p ^pV p . 

Definition 2.1. An Ehresmann connection A on P is a vertical-valued, one- 
form on P satisfying: 

j. A p : T p P —> V p is a linear map. 

2. A p is a vertical projection. That is, A p (v) = v Vv £ V p . 

■ 

The connection defines a horizontal space H p = ker A p at each point p £ P. The 
conditions in the definition imply T p P = V p ® H p and thus the connection gives a 
splitting of the tangent space at each point p into a vertical and a horizontal part 
(see again Figure 2.1). The union of these subspaces over all p is called the horizontal 
subbundle, i.e. H = U peP H p . 

Given an Ehresmann connection A, a point p £ P, and a tangent vector w £ 
Tn^B, define the horizontal lift of w to T p P as the unique tangent vector in H p that 
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projects to w under T p n. We call this lift hor p . The lift of w can be found by 

(2.1) hor p w = w — A p (w) 

where w £ T p P is an arbitrary vector satisfying T p tt(w) = w. 

Lemma 2.2. The map hor p is well-defined. 

Proof. Let ui, u 2 £ T p P be tangent vectors such that 

T p 7r(ui) = T p n(u 2 ) = U. 


Notice that 


T p tt (ui - u 2 ) = T p -k (ui) - T p -k (u 2 ) =u-u = 0 


and thus u\ — u 2 is a vertical vector. Let 


hor^ [«] = 2q — ApU\ and 


hor^ [u] =u 2 — A p u 2 . 


Then 


horp[u] - hor^ [u] = (ui - A p u{) - (u 2 - A p u 2 ) 
= (ui - « 2 ) - A p («! - tt 2 ) 

= (U! - tta) - (U! - u 2 ) = 0. 


□ 

It should be noted that while here we have defined the Ehresmann connection as 
a vertical-valued one-form and derived the horizontal space and horizontal lift, one 
can also begin with the definition of the horizontal space or the horizontal lift and 
define the other two objects. See [18] for details. 

2.2. Parallel transport and holonomy. Given an Ehresmann connection A 
on a fiber bundle, define parallel transport with respect to A along a curve lifted from 
the base space in the following way. Let b(t), t € [0,1], be a piecewise differentiable 
curve in B. The horizontal lift of b(t) with respect to A is the curve p(t) in P such 
that 7r= b(t) and that the tangent vector is horizontal for each t € [0,1]. 
We have the following proposition from [18]. 

Proposition 2.3. [18] Given a curve b(t), t € [0,1], in B and a point po € 
7T —1 (6(0)), there exists a unique locally defined horizontal liftp(t ) ofb(t) to P satisfying 
p(0) = po if P is a locally trivial fiber bundle. 

Proof. Since P is locally a trivial bundle we can write locally p = (b, /) and 
p = ( u , v) for b e B, f e F, u £ T^B, and v £ TfF. We can then write the connection 
one form as 


A(b,f)(u,v) = (0,v)+A(b,f)(u,0) 

= (0,u + A (b,f)u). 

(u,v) is horizontal if and only if A(b, f)(u,v) = 0 and thus for a horizontal tangent 
vector v = —A(6, f)u. If b(t) is a path in B denote p(t) = ( b(t),f(t )) where f(t) is the 
solution to the ordinary differential equation 


( 2 . 2 ) 


^(f) = -X(b(t),f(t))b(t), /(0) = / 0 where p 0 = (b 0 , / 0 ). 
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Then by local existence and uniqueness for ordinary differential equations this defines 
/(f) and thus p(t) for small t. If p(t) can be extended for all t £ [0,1] the connection 
is called complete. □ 

Given any curve b(t) in B, t € [0,1], and an Ehresmann connection A, the parallel 
transport operator Tb is defined as 

(2.3) T b : 7r —x (6(0)) -► tt 1 (^(1)), T b (p(0)) = p( 1) 

where p(0) £ 7r _1 6(0) and pit) is the horizontal lift of b(t ) with respect to A starting 
at p(0). 

By the uniqueness of the horizontal lift, is a bijection from 7 t _1 (6(0)) to 
7r _1 (6(l)) and by the smooth dependence of solutions of ordinary differential equations 
on initial conditions it is a diffeomorphism. 

Let bo be an arbitrary point of B and let C/ 0 be the set of all closed curves at bo, 
that is all b(t) such that 6(0) = 6(1) = 6 q. The diffeomorphism of 7r _1 (6o) onto itself 
given by parallel transport along 6(f) is called the holonomy of the path 6(f). Let 
<1%, be the collection of all parallel transport operators over Cb 0 and define the group 
operation as composition. 4>f, 0 then forms a group, called the holonomy group at 6o- 
(Assuming B is connected, it is easy to see that T(, 0 and T;,, are conjugate for any 
two 6o, 6i £ B. Thus if B is connected we have simply <3>, the holonomy group of the 
connection.) 

Definition 2.4. Given a bundle tt : P —> B, a connection on the bundle, and a 
closed curve b(t ) in the base space, the geometric phase is the holonomy along the 
curve b(t). m 

The geometric phase can be calculated by solving the ordinary differential equa¬ 
tion given in equation (2.2). 

3. The moving systems approach . Inspired by classical examples such as 
the Foucault pendulum and the ball in a hoop, one is led to consider the question 
of variations in phase space as a parameter in a mechanical system is slowly varied 
along a closed path in parameter space, as in a moving system. Examples in quantum 
physics, optics, and other settings [6, 26, 24] reveal that the essential calculation is a 
geometric one and is in fact captured by a holonomy. 

Marsden, Montgomery, and Ratiu have developed a modern geometric approach 
to understanding moving systems that utilizes the tools of Ehresmann connections 
on fiber bundles [18]. Here we present a brief synopsis of their approach. Let S 
be a Riemannian manifold and let M be the space of embeddings of a manifold Q 
into S. We think of S as the ambient space in which Q is being moved and of Q 
as the configuration space for a system of interest. A tangent vector to M at m 
is a map u m : Q —> TS such that u m (q) £ T m ^S. Given a tangent vector u m (q) 
one can construct a tangent vector to T q Q as follows. Relative to the metric on S, 
orthogonally project u m (q) to T rn ^m{Q) C ( T q m)(T q Q ), denote this vector u^q), 
and then pull-back u^(g) by Tm~ l to T q Q. This natural construction defines an 
Ehresmann connection on the product bundle tt : Q x M —> M as follows. 

Definition 3.1. [18] The Cartan connection omr : Q x M — > M is given by 
the vertical-valued one-form 7 C defined by 

(3.1) 7 c(q,m)(v q ,u m ) = (v q + (T~ 1 m o u^fq, 0). 


The Cartan connection induces a connection on p : x M -» M as follows. 
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Definition 3.2. [18] The induced Cartan connection on p :T*Q x M —> M 
is given by the vertical-valued one-form 7<, defined by 

(3.2) 7 oK ,m)(U aq ,u m ) = {U ag + X v{Um) {a q ), 0) 
where V(u m ) is the function defined by 

(3.3) (V(u m ))a q = a q -(T- 1 moul)(q) 

and is the Hamiltonian vector field ofV(u m ). m 

To separate the effects of the imposed motion on the system (as defined by the 
embeddings m t ) from the nominal dynamics (when the imposed motion is zero) we 
use the ideas of averaging. Abstractly, we assume we are given a left action of a Lie 
group G on T*Q. The average of the connection form 7 is defined by 

(3-4) <7) = P| J g 9*{l)dg 

where dg is a left Haar measure and |G| is the total volume of G. From this we have 
the following definition. 

Definition 3.3. [18] The Cartan-Hannay-Berry connection on p : T*Q x 

M —> M is given by the vertical-valued one-form 7 defined by 

(3.5) 7 K ,m)(U aq ,u m ) = (U aq +X {num)) (a q ), 0) 

where (•) denotes the average with respect to the action of the Lie group G. u In [18] 
Marsden, Montgomery, and Ratiu show that this is an Ehresmann connection. The 
horizontal lift of a vector field Z on M relative to 7 is 

(3.6) ( horZ)(a q ,m ) = (~X {nz(m))) (a q ), Z(m)) . 

Definition 3.4. [18] The holonomy of the Cartan-Hannay-Berry connection is 
called the Hannay-Berry phase for a moving system. ■ 

4. The rotating, vibrating ring . In this section we derive the Hannay-Berry 
phase for the vibrating ring. Using the moving systems approach we find an explicit 
formula for the phase shift under linearizing assumptions and show that this result 
matches that of G.H. Bryan [8]. In this work we are interested in the effects of the 
imposed rotatory motion and as a consequence choose to simplify the analysis of the 
ring dynamics by assuming the ring has no cross-sectional area. This choice also allows 
a direct comparison to the results derived by Bryan. A more comprehensive treatment 
based on the geometrically exact theory of rods could be developed to understand the 
detailed dynamics of the ring itself (see, for example, [20]). 

Consider a thin ring of length L and line density a. The body is given by B = 
{b : b £ [0, L}}. Let 9 be the mapping given by 

9 : B -► S 1 



allowing us to parametrize the ring by 9 e [0,27r]. We define the reference configu¬ 
ration to be a circular ring of radius a centered on an inertial reference frame and 
(w(9),j(9)) to be the radial and angular deformations from this reference respectively. 



THE BERRY-HANNAY PHASE OF THE RING GYROSCOPE 


To maintain integrity of the ring we require «,'(0) = w(2tt) and 7 ( 0 ) = 7 ( 277 ). In stan¬ 
dard cylindrical coordinates the configuration of the ring is given by (a+w(8),6+'y(0)). 
Let C be the space of all smooth deformations of the ring. (We do not discuss here 
the explicit infinite dimensional manifold structure for C and associated structures, 
although it is standard as in [17]). Since we are interested in imposed rotational 
movements of the ring (as a sensor), we split j(8) = ip + a{6) with a(0) = a(27r) 
where ip, independent of 6, is a global rotation. 

We now use the following argument of Rayleigh [23]. Since the ring is thin the 
forces resisting bending are small in comparison to those which resist extension. In 
the limiting case of an infinitely thin ring the flexural vibrations become independent 
of any extension of the circumference as a whole and one may assume that each part 
of the circumference retains its natural length throughout the motion. Under this 
condition we say the ring is inextensible. Viewing the deformed ring as a curve in 1R 2 , 
a point on the curve is given in Cartesian coordinates by 


(4.1) 


( x(9) \ ( (a+ «;(*)) cos(0.+ 7(0)) \ 
V y(°) ) V (a + ■*#)) sin(0 + 7 ( 0 )) ) 


Equating the lengths of an arbitrary section of the circumference of the reference 
configuration to the length of the same section in the deformed configuration yields 


L' aM= l'\ 

/ / dxP 

2 / 0 \ 2 


P02 

(4 ' 2) -I, \ 

j (a + i 

'< s » 2 ( i+ §) 

2 ( dw \\o 
+ de 

Since and #2 are arbitrary we 

i have 



O = (o + U, ) 2 (i + I 0 + 

fdw\ 

\de) 

2 

— a 2 



From here on we assume the deformations are small and so we keep only terms to 
first order in equation (4.3). The inextensibility condition then requires that 



From the above the space C is given by 
(4.5) C = {(ip,a)\ip £ S' 1 , a : S 1 —> S' 1 ,a(0) = a(27r),a smooth}. 


Any W £ has the form (with overdot denoting the partial derivative with 

respect to time) 


<«> w =(m)- 

The kinetic energy is easily verified to be 
(4.7) KE=\ jT [(l - (> + 2 iPd(6) + a\6)) + (^j j oa 3 d0 
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where equation (4.4) has been used to express w in terms of a. This defines an inner 
product on C given by 


(Wi 


,w 2 ) = J (v>iV >2 +^ 2 ( 21 ( 0 ) + ipid 2 (9) +ai(9)a 2 {0)^ 


(4.8) 


ddi dd 2 

f ~W~dO 


As in Bryan, we take the potential energy due to the bending of the ring to be 
proportional to the change in curvature squared of the ring. That is 

(4.9) V-fjf (Ka(9) - n a =o(0)f ade 

where K a (9) is the curvature of the ring at the material point 6 under the deformation 
a and (3 is a material constant. For a curve (r(t), 4>{t)) defined in polar coordinates, 
the curvature is given by 


(4.10) 


_ (2 r 2 4> + rr<j> - rr<j> + r 2 <p 3 ) 
(f 2 — r 2 o 2 ) 5 


In the 9— parametrization, the configuration of the ring under the deformation 
(a, w ) is given by the curve 


t 411 ) r ( 6 ') - a, + w(9), 

<j>(9) = 9 + ip + a(9). 

Using equation (4.11) in (4.10) to express the curvature of the ring under the de¬ 
formation ( w,a ), we find (keeping only terms to first-order in the numerator and 
denominator) 


(4.12) 


( — + o? + 2 aw + 3a 2 ^) 
(a 2 + 3 a 2 w + 3a 3 |f) 


Writing the curvature in terms of a alone using the inextensibility condition of equa¬ 
tion (4.4) yields 


(4.13) 



From this we have K a =o(9) = i. Using these results in the potential energy of equation 
(4.9) yields 


(4.14) 


V = 


p_ r 2n ra 3 ^ 

2a J 0 [W + 


da' 2 
89_ 


dJ9. 


In Bryan’s original work the potential energy included also a term capturing the 
work done in stretching the ring and the work done against an attracting force which 
he introduced to separate the effects of the centrifugal force from the remaining terms. 
These two terms are related by a simple equation involving the rate of the imposed 
rotation and at the conclusion of his analysis Bryan chooses the attracting force so as 
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to cancel the tension, leaving only the work done in bending. Here we take a simpler 
approach, similar to Rayleigh, and omit those terms at the outset. 

The standard Lagrangian function is defined to be the kinetic minus potential 
energies and is given here by 

L = \[' [(' - wf + 2 + “ a < e >) + (§)] 

Consider now the following action <E> S of S ' 1 on C. 

(4.16) <& g (ip,a) = (ip + g,a). 

We have the following lemma. 

Lemma 4.1. (C, (•, •), V, S 1 ) is a simple mechanical system with symmetry where 
the action of S 1 onC is given by equation (4.16). (For a definition and discussion of 
simple mechanical systems with symmetry see Section 4.5 of [1] J 

Proof Immediate since both (•,•) and V are invariant under the given action of 
S 1 on C. □ 

Since the given action is both free and proper, the reduced space Q = C/S 1 , given 
by 

(4.17) Q = {a : S 1 —» 5 1 |a(0) = a(27r),a smooth} 


is also a manifold. To fix notation in relation to the general theory presented in 
Section 3, we note that Q = C/S 1 is the configuration space for the ring and S' = C is 
the ambient space in which Q is moved. To slowly rotate the ring we set ip = ipo + Clt 
(identifying ip = 0 with ip = 2 tt) for some small ft and some fixed initial offset ipo so 
that the embedding from Q to S is given by 

(4.18) m t (a(0)) = (ip 0 + Qt,a(9)). 


4.1. The nominal dynamics. The nominal dynamics is given by setting fl = 0 
in equation (4.18). Applying this to equation (4.15) yields the nominal Lagrangian. 


—r {?[(-£ 


V^.fSV 


e) ^\oe) 

]> 

The action integral for this Lagrangian is defined to be 


(4.19) 


13 rd 3 a dal I J/t 

2a [oe 3 + do\ \ de 


(4.20) 


m (a, a)dt. 


J(a, a) = J . 

as for this 
.g. [ 1 ]) w 

f b 

6J(a,a) = 6 L Nom (a,a)dt = 0 


The Euler-Lagrange equations for this system are found by applying Hamilton’s 
principle of critical action (see, e.g. [ 1 ]) which states that 


(4.21) 
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for all variations among paths r)(t) in Q with fixed end-points. Applying the variation 
yields 

SJ(a,a) = — J(a + erj,a + erj) 

dc | e _Q 



Using integration by parts repeatably on the space variable and the fact that for any 
element £(■) e Q we have £(0 ) = £(27t), the variation of the action can be rewritten 
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convenient set of coordinates, namely the coeffecients of a, a expressed in a Fourier 
basis. We first express the nominal dynamics in these coordinates and in the following 
section do the holonomy calculations in the same coordinates (after truncation of the 
Fourier series). In these coordinates a has the form 


(4.24) a(0) = 22 \ A k cos (kd) + B k sin(fcfl)]. 

k=i 

The deformation a(0) is not allowed to contain any global rotations and so the con¬ 
stant coefficient is set to 0. Inserting this expression for a into the equation of motion 
(4.23) results in the equation 


0 = '22 |cra 3 (l + k 2 ) ^A k cos (k9) + B k sin(fc0)J 

(4.25) + 1(k 6 - 2/e 4 + k 2 ) [A k cos (k6) + B k sin(fc0)] J . 

Collecting terms in cos(0) and sin(d) and setting them separately to zero gives the 
following set of ordinary differential equations for the Fourier coefficients. 


(4.26) 

(4.27) 


A k 

B k 


P k 2 (k 2 — l) 2 „ a 
era 4 k 2 +1 k ~ 

P k 2 (k 2 - l) 2 A 
era 4 k 2 +1 k 


-Vk A k, 


-vl B k 


which defines the frequencies r) k . This result is in agreement with a derivation of 
Rayleigh [23] and defines for each k a pair of uncoupled oscillators with common 
frequency r] k . The solution to this system is given by 


(4.28) 

(4.29) 


A k (t) = A k cos (ij k t) + — sin(%f), 
Vk 

B k (t) = B k cos(r?fcf) + — sin (r] k i) 
Vk 


where A k , A k , B k , and B k are given by initial conditions. The Hannay-Berry phase 
is defined as the holonomy on a trivial bundle involving the cotangent space of the 
system. It will prove useful, then, to have the time evolution of the conjugate momenta 
for the nominal system. By inserting the Fourier expansion for a into the nominal 
Lagrangian, equation (4.15), and applying the Legendre transform we obtain 


(4.30) 

(4.31) 


PAk = ^T = ^ + 

dA k 

PB k = • = (1 + k 2 ) 

dB k 


i 3 nA k , 


Thus the solution to the nominal system expressed on the cotangent bundle is given 
by 


(4.32) 


A k (t) = A k cos(rj k t) + 


(1 + k 2 )aa 3 nr] k 


sin (r?fcf), 
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(4.33) 

PA k (t ) = -(1- 

1- k 2 )aa 

3 nr] k A k sin (<q k t) + p Ak cos (%t,) 

(4.34) 

B k (t) = B k cc 

»s {r]kt) H 

1- /I , l2 s k o sm (r] k t), 

(1 + k 2 )aa 3 Trr] k 

(4.35) 

PB k ( t ) = -(1 - 

1- k 2 )aa 

3 nT] k B k sin (rj k t) + p Bk cos (r) k t). 


4.2. The Hannay-Berry phase. From equation (4.18), the velocity vector of 
the motion in S is 

(4-36) |(m t (a(0))) = (O,a(0)) + (f2,O). 

From this we see that the vector field corresponding to the imposed motion is 

(4.37) Z(m t (a(6))) = ( q ) ' 


For ease of notation define Z = Z t {rn t {a{0))). The projection of this tangent vector 
to T mt ( q yrn t (Q) with respect to the metric of S is given by Z T = Z — Z 1 - where 
(^- L ,X) = 0 \/X e T mt ^m t {Q). From equation (4.6), any vector X e T mt ^m t (Q) 
has the form 

(4.38) x= (y 

where Y e T q Q. Then, from equation (4.37) and the fact that Z T e T rrit ( q) rn t (Q), we 
can write Z T = Z — Z 1 - as 


(4.39) 



for some Y z t,Z ±2 6 T q (Q) and Z 3 - 1 e T^S 1 . From equation (4.39), Z- 1 - 1 = fi. 
Applying equation (4.8), the orthogonality condition states 

o = ((fi,£ ± 2 ),(o,y)) 

(4.40) 

In what follows, we express the orthogonality condition of equation (4.40) in the 
Fourier basis and thus derive an explicit formula for Z ±2 (see (4.51) and (4.52) below). 
Using the Fourier series representation for a, a tangent vector Y € T q Q has the form 


-£) 


(4.41) y = l y A k cos(fcfl) + Y Bk sin(fc6>)] 

Adopting the notation c (0) = cos (0) and s (6) = sin(d), the orthogonality condition is 
given by 
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Using the follow 

ing identities 

(4.43) 

J * c(k6)c(W)dO = 

( 4 - 44 ) J 

■27r 

c(k0)c(W)c{m0)c{nO)d0 = 

(4.45) J 

r 2n 

s(k0)s(l0)s(m0)s(n0)d0 = 

(4.46) J 

-2tt 

c(ke)c(W)s(me)s{n0)d0 = 

where 


(4.47) 5^ = 

r i if fc = f 
\ 0 otherwise 

and the fact that all other combinations of 
integrate to 0 over [ 0 , 2ir\, the orthogonality 


'E [( z t +t Alz *i + 

+ k 2 Z)A Y Ak 


This holds for every Y 6 T q Q and so for every k we have 

(4.49) Zjg + k *Alzil + ^B\Z\1 - ^ A k B k Z 4f - 2f2fcB fc + k 2 Z^* = 0 

01.2 1.2 7,2 

(4.50) + 4 A\Z%. + -B 2 Zt h - 2 A k B k Z% + 2SlkA k + k 2 Z^ = 0. 

Solving these coupled equations for Z A * and Z^ yields 

„ C1 , ^t 2 0 \m + k 2 + s fA 2 + * a B 2 )B k - k 3 A 2 B k ] 

(451) ^ ° [ D k (A k ,B k ) J ’ 

<-> z ,- n [ ^ + p + ¥4^-^ i 


(4.53) 
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and so, representing the tangent vector by its coefficients at each k, 


q 2fc(l +k 2 +^Al+^Bl)A k -k 3 A k Bl 


Inserting equation (4.54) into equation (4.39) gives the projection of the tangent vector 
of the imposed motion onto T mt ^m t (Q) to be 


(4.55) = 


!fc(l+fc 2 + VA^ + 5j-Rg)A fc -fc 3 A fc Bg 


The pull-back of Z T to T q Q by [ Tm] 1 is given by 


2fe(l+fc 2 + ^f— B^)Bk~ fe 3 A^Bk 


(4.5 6)Z(q) = [Tm\- 1 Z T = 


where Z(q ) is defined for ease of notation. Recalling that the deformations are 
assumed to be small, the above expression is expanded in a Taylor series about 
Ak = Bj- = 0 V k and only the first order terms kept. This yields 


Let qk = ( Ak,Bk ) so that the coordinates on Q are q = {g&ffTr The conjugate 
momenta can then be written pk = ( PA k ,PB k )• Define 

/. ro'i 7 , , A 2kn (o -1 \ f Ak \ 2fcn 

(4.58) = o ) ( B J = TTP S(1)ft 

which also defines the skew symmetric matrix 5(1). With this definition the projected 
vector can be expressed as Z(q) = {Zk{Qk)\ 12a- To avoid technical difficulties we 
assume that only N Fourier modes are active where N is some positive, finite number. 
Using this assumption, the function V(Z) defining the horizontal lift relative to the 
induced Cartan connection (as in equation (3.3)) is given by 
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Define Q k = {(Ak,Bk) € IR 2 } so that Q = U k =iQk- The configuration space is 
then the Cartesian product of N copies of IR 2 and each coordinate Qk and conjugate 
momenta Pk can be identified with a vector in ]R 2 . Extend these vectors to IR 3 by 
letting the third coordinate of each be zero. Let x 3k be a unit vector at the origin of 
the k th copy of IR 2 along this third direction. With these identifications, in equation 

(4.59) we replace pk ■ S(l)qk by (qk x p k ) ■ x 3k . Define 

(4.60) 4 = ( q k x p k ) • x 3k 
so that 

J5L 2k 

( 4 - 61 ) P{Z)(q,p) = ^ 

Let T be the subset of C°° functions on T*Q defined by 

(4.62) T = ^f(q,p)\f(q,p) = ^ a k f k (A k , B k ,PA k ,PB k ), a k € ]R,/ fe STOOotft.|. 

Since the time solution to each Fourier mode for the nominal system is periodic, we 
can define an average on T with respect to the flow by 

N 

(4-63) < f)=J2a k (f k ) k 


where 

(4-64) (fk)k=^J o Vk {^Tfkdt 

where $ is the flow corresponding to the k th Fourier mode of the nominal system. 
From equations (4.60) and (4.61) it is clear that P(Z) € T. We have the following 
useful lemma. 

Lemma 4.2. 4 constant along the trajectories of the nominal system. 

Proof. 

dh d 

n = *(«*»>•*«, 

= j t (PB k A k -PA k B k ) 

= i>B k A k + PB k Ak - i>A k B k ~PA k B k 

= (1 + k 2 )aa 3 Trr]k ( B k A k + B k A k - A k B k - A k B 
= (1 + k 2 )aa 3 Trr]k ( B k A k - A k Bkj 

where in the second to last step we have used the definition of the conjugate momenta 
for the nominal system in equations (4.30,4.31). From equations (4.26,4.27), along 
the trajectories of the nominal system we have 

( B k A k - A k B k ) = -vl (A k B k - A k Bk) = 0 . 


(4.65) 
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□ 

The average of V(Z) over the nominal dynamics is then 
N ?k E 2k 

(4.66) = •£ = E r Vk^ h - 

k=1 k=1 

Noticing that this function depends only the coordinates ( qk,Pk ) through /; c , we move 
to the coordinates for the averaged dynamics in phase space defined by (cpk, Ik, Pk,Pp k ) 
where pk = ((A| + Here (pk is conjugate to Ik and p Pk is conjugate to pk- 

With these coordinates the horizontal lift of ft relative to the Cartan-Hannay- 
Berry connection (as in equation (3.6)) is given by 

(4.67) t ,n) - 

The geometric phase is neatly split into a phase change in each (pk independently. If 
the loop in M is parametrized by t € [0, T] where T = is the time to complete one 
full revolution of the ring, then 

(4.68) ^ ndt = -TT^ QT = - 2 "'[r^fcs] - 

After one full rotation of the ring each vector qk has been rotated by the angle A (pk- 
In practice one expects to use only one mode under a resonant drive (as in [21, 22]). 

Note that we perform an average with respect to the flow of the nominal dynamics 
on a special class of functions as in equation (4.62). This agrees with a group (here 
S 1 ) average as in Section 3 (Definition 3.3) when we restrict the nominal dynamics 
to a single mode. 

5. A comparison with the results of Bryan . In [8] Bryan uses classical 
variational techniques to derive the equations of motion for a thin ring of radius a 
undergoing a steady rotation about its central axis with angular velocity f 1. His 
analysis uses two polar coordinate systems, the first fixed in space and the second 
rotating with angular velocity f2. If in the undeformed state the coordinate systems 
are given by (a, (p) and (a, 6*) we have 

(5.1) 4>=0-(It 

and 6 is constant for any particle of the ring. Let the tangential and radial displace¬ 
ments of a particle of the ring be given by v and w respectively so that the new 
polar coordinates are (a + w,(j>+ v/a) and (a + w,9 + v/a) in the two systems. The 
deformations v, w are assumed to be small. The assumption of inextensibility yields 


as before. As discussed in the previous section, Bryan includes work done against 
the tension, T, to stretch the ring and against an attractive force, //. To match 
his derivation with the model we have chosen, we set these terms to zero. Taking 
variations on the total energy and setting them to zero, we find the following equation 
of motion 
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Bryan then assumes that the deformations are of the forir 


Inserting this into the equation of motion yields the following two solutions. 


_2 _ /? k 2 (k 2 - l) 2 _ fl 2 k 2 (k 2 - 3) 

° Jk era 4 1 + k 2 (1 + k 2 ) 2 

( 57 ) „2 » 2 t 2 (* 2 - 3 ) 

l 5 ') It (1 + 02)2 • 

Notice that by retaining the terms in fl 2 there is a slight decrease in the frequency 
of vibration from %. This can be understood as a “softening” of the material and 
corresponds to the spurious softening that occurs in the theory of rotating rods if 
the models are linearized prematurely, that is before the effects of external rotation 
are considered. The geometrically exact theory handles this issue properly and we 
believe that extending the model of the ring using this theory will prove useful. (We 
note also that if the attractive force /x introduced by Bryan is kept in the equations 
of motion and set equal to fl 2 then the resulting system shows an increase in the 
frequency of vibration due to the imposed rotation. See [13] for comments on similar 
ad hoc methods in the theory of rotating beams.) 

Assuming the amplitude A in the two solutions is the same, the final solution in 
the fixed frame is given by 


Bryan then recognizes that this corresponds to an oscillation with 2 k nodes where the 
position of the nodes precess in retrograde around the ring with angular velocity 


e this in the rotating system we have 

(pTT- 1 ) 0 —PTI n 

one rotation the nodes have precessed in the moving frame by 


To compare this to 
ingle mode so that 


results in Section 4 we must first restrict our solution t 


a(6,t) = Ak(t) cos (k0) + Bk(t) sin (kO). 
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When ( Ak,Bk ) is viewed as a vector in 1R 2 , the effect of the geometric phase is seen 
to be a rotation of this vector about the origin where the counter-clockwise direction 
is defined to be positive. Using equation (4.68), the rotated vector at the end of one 
revolution of the ring is given by 


( MT) 
V MT) 


(5.13) 


cos (-27^) -sin (-27^) \ / ^ fc(0 ) \ 

sin (-277^) cos (—27r r |^) J V M°) ) 

^fc(O) cos (-2?TiHr) - B k (0 ) sin \ 

^4fc(0) sin 27rJ + B k {0 ) cos J 


Inserting this into equation (5.12) and simplifying we get 


(5.14) 


a{9,T) = A fc (0)cos^ 
+B k ( 0) sin 


yh^TTp]) 

(‘h^ITp]) 


which is course the same solution with the nodes rotated by — 2n J, agreeing 

with Bryan. 

6. Nonlinear corrections . We now turn to an investigation of corrections to 
the geometric phase based on the nonlinear terms in the vector field Z(q) given in 
equation (4.56). It is worth noting that these arise due to the configuration-dependent 
quadratic form defining the kinetic energy. We proceed by keeping higher-order terms 
in the Taylor expansion of Z(q). The second-order terms in this expansion can be 
shown to be zero. To third order the vector field is 


( 6 . 1 ) 


/ 

Z(q) = 

V 


fc 2 (2 +k)A 2 k B k 9fc 3 .Bgl'l Ar 


Mm* 


k 2 {2+k)A k B 

uw 


\ 


) 


Define the matrix U(q k ) by 


( 6 . 2 ) 


U{q k ) 


With this, equation (6.1) can be written as 

z(q) = {tTP*' 1 ’® - + 

qr.3Q 1 N 

( 6 - 3 ) -(IW^^^Lr 

and the function defining the Hannay-Berry phase is given by 

nm.P )=E • s(D« - ■ v(kW)<* 
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(6.4) 



2kCl 
1 + fc 2 


4 


% + ' S ^ U ^ Qk 


k 2 (2 + k)n 

(1 + k 2 ) 2 


Pk ■ U(q k )S(l)q k 


Qb 3 fl 

~ (1 + k 2^ P k - S ( 1 M qk )9 k 


where in the second step the definition of I k from equation (4.60) has been used. 
To determine the Hannay-Berry phase we need to find the average of equation (6.4) 
over the nominal dynamics. In Section 4 we have shown that the first term in square 
brackets in the sum in equation (6.4) is a constant along trajectories of the nominal 
system. The second and third terms, however, are not constant and their averages 
need to be explicitly calculated. For the second term the average of the k th element 
of the sum is given by 

(Pk • U(q k )S(l)q k ) k = | k -j Vk [- p Ak (t)Al(t)B k (t ) 

(6.5) +p Bk ( t ) A k (t) B\ (t)] dt. 


Using the solution to the nominal system given in equations (4.32-4.35) and making 
a change of variables in the integration, this becomes 


{Pk ■ U(q k )S(l)q k ) k 


1 f 2n 

— J {[(1 + k 2 )aa 3 Trr) k A k sin(t) - p Ak cos(t)] 


( 6 . 6 ) 


• cos (t) + 

+ [-(1 + k 2 ) 

■ \ A k cos( t) 4 


27T 

{l + k 2 )aa 3 TTi lk \ [ 
Ta 3 TTT] k B k sin(f) + p Bk cos(f)] 
PA k sin(t) If 


(1 + k 2 )aa 3 nr]k\ 


PB k sin(f) 1 
(1 + k 2 )aa 3 irq k \ 

PB k 

(1 + k 2 )aa 3 nri k 


where, through a standard abuse of notation in averaging, the hats have been dropped 
on the initial conditions. The integration identities in equations (4.43 - 4.46) can be 
used to write the above expression as 

{Pk ■ U(q k )S(l)q k ) k = i [(Alp Bk - B%p Ak - A 2 k B k p Ak + A k B k p Bk ) 

(A k p 2 Ak p Bk ~ B k p Ak p% - B k p\ + A k p% \ 1 

(8 - 7) + (- ( (i + F)^) a - )\- 

Now 


A kPB k ~ B kPAk ~ A l B kPA k + A k B k p Bk = ( A kPB k ~ B kPA k )( A l + B l) 
(6.8) =h{A 2 +B 2 k ) 

and 


A kPA k PB k - BkPA k P% k - B kP Ak + A kPBk = ( A kPBk ~ B kPA k )(p Ak + P %J 

= h(p 2 Ak +pl k ) 


(6.9) 
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and so equation (6.7) takes the form 


(6.10) (p k • U(q k )S(l)q k ) k = ^ I A 2 + B 2 + 


PA k +Pb , 


((1 + A: 2 )cra : 


^_ ] 

•nrik) 2 ) 


Following the same procedure, the average of the k th element in the sum of the 
third term is 

{p k ■ S(l)U(q k )q k )k = | ^ ( PB k (t)A 3 k (t ) - p Ak (t)Bl(t))dt 

= g [( A k B k p Bk - A 2 k B k p Ak + A\p Bk - B\p Ak ) 

A kP% k - B kP Ak + A k p 2 Ak PB k - B k p Ak p% k 

+ ((1 + k 2 )aa 3 nr]k) 2 


( 6 . 11 ) 


A 2 + B 2 + 7 


( 6 . 12 ) 


((1 + k 2 )aa 3 nr] k ) 2 
Using equations (6.10,6.11), the average of V(Z) is given by 

jl r 2 k 

k 2 + 14A; 3 ( A l + B k+ ((l+fe^cr^V) 2 


2(1 +A: 2 ) 2 


From equations (4.32,4.34) we see that the term in parentheses is the average 
of (A‘ 2 k + B' k ) over the nominal dynamics. We move to the averaged coordinates 
(< p k ,I k ,p k ,P Pk ) as in the comments following equation (4.66). In these coordinates, 
the average of V(Z) has the form 


(6.13) 


(V(Z)((q,p) = Y, ni k YT 


k 2 + 14fc 3 
2(1 + k 2 ) 2f> 


The lift to third-order of ft with respect to the Cartan-Hannay-Berry connection 
is given by 


(6.14) 


' 2k 

k 2 + 14fc 3 

1 + k 2 

2(1 +A: 2 ) 2 ' 

-fll k pk 

Ik 2 + 14A 3 ! 

[(l + k 2 ) 2 \ , 


' d 

d<t> k 


From equation (6.14) we see that both I k and p k are constant. Thus 
a r T n \ ™ k 2 + Uk 3 2 ] 

afc=_ X a [m?-w^W’ k \ dt 

k 2 + 14k 3 


(6.15) 


[l + k 2 2(1 +A; 2 ) 2 


Pk 
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Notice that the third-order terms act to reduce the rate of nodal rotation and thus 
the sensitivity of a vibrating ring gyroscope cannot be increased by increasing the 
amplitude of vibration and using the nonlinear effects. 

In contrast to the earlier calculation where we kept in Z only the terms linear in 
configuration variables, in the present nonlinear setting the imposed rotation causes 
not only a precession of the nodes of vibration but also a drift in the momentum 
conjugate to pk■ In practical devices the ring is driven into a single mode of oscillation 
and the imposed rotation sensed by measuring the drift rate of the nodal points of the 
vibration. Thus the effect of the second term in equation (6.14) will be compensated 
for by the drive electronics. 

It is interesting to ask how large the nonlinear effect on the drift rate of the nodal 
points of the vibrations is in a typical device. The micromachined ring gyroscope 
of Putty and Najafi [22] utilizes a ring of radius a = 500 pm placed into elliptical 
vibration so that k = 2 with a radial deformation amplitude of 0.15 pm. From 
equations (4.4,4.24), the radial deformation for this ring is 

(6.16) w(9) = 2 a [A 2 sin (2 6) - B 2 cos(20)]. 

Let t = 0 be the time at which the maximum radial deformation is attained and 
let 9 = 0 be the location on the ring of the maximum radial deformation at time 
t = 0. From these definitions we have the initial conditions A 2 = pa 2 = Pb 2 =0 and 
B‘2 = — 2 ( 500 ■ Inserting these values into equation (6.15) yields 

(6.17) A </>2 = -27r - 2.61 x 10“ 8 j . 

For the normal operation of this device, then, the nonlinear effects are seven orders 
of magnitude smaller than the first-order terms. 

7. Conclusions and future work . In this paper we have applied the moving 
systems approach developed by Marsden, Montgomery, and Ratiu to the rotating, 
vibrating ring and showed that the precession of the nodes with respect to the ring can 
be understood as the Hannay-Berry phase. We have also showed that under the same 
linearizing assumptions, our result matches that of Bryan’s original analysis. Using 
the inherently nonlinear nature of the moving systems approach we then calculated 
the effect of the imposed rotation on the nodal precession to third-order. These higher- 
order terms serve to reduce the sensitivity of the device and we therefore conclude 
that the best performance of these devices is obtained by operating them in the linear 
regime. 

There is an underlying assumption that the rate of rotation of the ring is slow 
with respect to the frequency of vibration. Under this assumption the character of the 
nominal dynamics is unaffected by the imposed motion (in the linearized setting) and 
the effect of the rotation is completely captured by the nodal precession. At rotation 
rates that are not small with respect to the rate of vibration the dynamics are severely 
affected and these techniques do not apply. However even at slow rotation rates the 
effect on the vibration is non-zero. In a future paper we will report on recent work in 
which we extend the theory of moving systems using Hamiltonian normal form theory 
to account for the non-adiabatic nature of the imposed motion (see also [ 2 ]). 
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